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Within the framework of the kinetic energy driven superconductivity, the electronic structure 
of bilayer cuprate superconductors in the superconducting state is studied. It is shown that the 
electron spectrum of bilayer cuprate superconductors is split into the bonding and antibonding 
components by the bilayer splitting, then the observed peak-dip-hump structure around the 0] 
point is mainly caused by this bilayer splitting, with the superconducting peak being related to 
the antibonding component, and the hump being formed by the bonding component. The spectral 
weight increases with increasing the doping concentration. In analogy to the normal state case, both 
electron antibonding peak and bonding hump have the weak dispersions around the [n,0] point. 
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I. INTRODUCTION 

The parent compounds of cuprate superconductors are 
the Mott insulators with an antiferromagnetic (AF) long- 
range order (AFLRO), then via the charge carrier doping, 
one can drive these materials through a metal-insulating 
transition and enter the superconducting (SC) domei££. 
It has become clear in the past twenty years that cuprate 
superconductors are among the most complex systems 
studied in condensed matter physics^^. The compli- 
cations arise mainly from (1) a layered crystal structure 
with one or more CuC>2 planes per unit cell separated by 
insulating layers which leads to a quasi-two-dimensional 
electronic structure, and (2) extreme sensitivity of the 
physical properties to the compositions (stoichiometry) 
which control the carrier density in the CuC>2 plane^^. 
As a consequence, both experimental investigation and 
theoretical understanding are extremely difficult. 

By virtue of systematic studies using the angle-resolved 
photoemission spectroscopy (ARPES), the low-energy 
electronic structure of cuprate superconductors in the 
SC state is well-established by now 2 -^, where an agree- 
ment has emerged that the electronic quasiparticle-likc 
excitations are well defined, and are the entities par- 
ticipating in the SC pairing. In particular, the lowest 
energy states are located at the [tt, 0] point of the Bril- 
louin zone, where the d-wave SC gap function is max- 
imal, then the most contributions of the electron spec- 
tral function come from the [tt, 0] points. Moreover, 
some ARPES experimental results unambiguously estab- 
lished the Bogoliubov-quasiparticle nature of the sharp 
SC quasiparticle peak near the [tt, 0] points, then the 
SC coherence of the quasiparticle peak is described by 
the simple Bardeen-Cooper-Schrieffer (BCS) formalism 6 . 
However, there are numerous anomalies for different fam- 
ilies of cuprate superconductors, which complicate the 
physical properties of the electronic structure 2 -^. Among 
these anomalies is the dramatic change in the spectral 
lineshape around the [tt, 0] point first observed on the bi- 
layer cuprate superconductor Bi2Sr2CaCu20s+5, where 
a sharp quasiparticle peak develops at the lowest bind- 
ing energy, followed by a dip and a hump, giving rise to 
the so-called peak-dip-hump (PDH) structure in the elec- 



tron spectrum^^. Later, this PDH structure was also 
found in YBa2Cu307_^ and in Bi 2 Sr2Ca2Cu30io+^i. 
Furthermore, although the sharp quasiparticle peaks are 
identified in the SC state along the entire Fermi surface, 
the PDH structure is most strongly developed around the 
[tt,0] pointiLSAiaJl. 

The appearance of the PDH structure in bilayer 
cuprate superconductors in the SC state is the mostly 
remarkable effect, however, its full understanding is still 
a challenging issue. The earlier work a 2 : 12 gave the main 
impetus for a phenomenological description of the single- 
particle excitations in terms of an interaction between 
quasiparticles and collective modes, which is of fun- 
damental relevance to the nature of superconductivity 
and the pairing mechanism in cuprate superconductors. 
However, the different interpretive scenario has been 
proposed 2 * 1 ^. This followed from the observation of the 
bilayer splitting (BS) for both normal and SC states 
in a wide doping rang o 14 i 15 ' 16 . This BS of the Cu02 
plane derives the electronic structure in the bonding and 
antibonding bands due to the present of Cu02 bilayer 
blocks in the unit cell of bilayer cuprate superconduc- 
tors, then the main features of the PDH structure is 
caused by the BS 13 i 14 ' 15 i 16 , with the peak and hump 
corresponding to the antibonding and bonding bands, 
respectively. Furthermore, some ARPES experimental 
data measured above and below the SC transition tem- 
perature show that this PDH structure is totally unre- 
lated to superconductivity^ 4 .. The recent ARPES exper- 
imental results reported by several groups support this 
scenario, and most convincingly suggested that the PDH 
structure originates from the BS at any doping levels 1 ™ 
To the best of our knowledge, this PDH structure in bi- 
layer cuprate superconductors has not been treated start- 
ing from a microscopic SC theory. 

Within the single layer t-t'-J model, the electronic 
structure of the single layer cuprate superconductors in 
the SC state has been discussed 1 ^ based on the frame- 
work of the kinetic energy driven superconductivity-^, 
and the main features of the ARPES experiments on 
the single layer cuprate superconductors have been repro- 
duced, including the doping and temperature dependence 
of the electron spectrum and quasiparticle dispersion. In 
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this paper, we study the electronic structure of bilayer 
cuprate superconductors in the SC state along with this 
line. Within the kinetic energy driven SC mechanism 19 , 
we employed the t-t'-J model by considering the bilayer 
interaction, and then show explicitly that the BS occurs 
due to this bilayer interaction. In this case, the elec- 
tron spectrum is split into the bonding and antibond- 
ing components by this BS, then the SC peak is closely 
related to the antibonding component, while the hump 
is mainly formed by the bonding component. In other 
words, the well pronounced PDH structure in the electron 
spectrum of bilayer cuprate superconductors is mainly 
caused by the BS. Furthermore, the spectral weight in 
the [it, 0] point increases with increasing the doping con- 
centration. In analogy to the normal-state cas o 14 ' 20 ' 21 ' 22 , 
both electron antibonding peak and bonding hump have 
the weak dispersions around the [tt,0] point, in qualita- 
tive agreement with the experimental observation on bi- 
layer cuprate superconductors in the SC state 2 -^ ^ 10 ' 11 . 

The paper is organized as follows. The basic formal- 
ism is presented in Sec. II, where we generalize the 
kinetic energy driven superconductivity from the previ- 
ous single layer cas o 18 i 19 to the bilayer case, and then 
evaluate explicitly the longitudinal and transverse com- 
ponents of the electron normal and anomalous Green's 
functions (hence the bonding and antibonding electron 
spectral functions). Within this theoretical framework, 
we discuss the electronic structure of bilayer cuprate su- 
perconductors in the SC state in Sec. III. It is shown that 
the striking PDH structure in bilayer cuprate supercon- 
ductors is closely related to the BS. Finally, we give a 
summary and discussions in Sec. IV. 



II. FORMALISM 

It has been shown from the ARPES experiments 2 -^ 
that the two-dimensional t-t'-J model is of particular 
relevance to the low energy features of cuprate super- 
conductors. For discussions of the physical properties of 
bilayer cuprate superconductors, the t-t'-J model can be 
expressed by including the bilayer interactions as, 
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supplemented by an important on-site local constraint 
J2 a Cj aa Ci a(T < 1 to avoid the double occupancy, where 
f) = ±x, ±y representing the nearest neighbors of a given 
site i, t = ±x± y representing the next nearest neighbors 
of a given site i, a = 1,2 is plane index, C\ aa (Ci aa ) 
is the electron creation (annihilation) operator, Si a = 
Cl a crCi a /2 is the spin operator with the Pauli matrices 
er = (cr Xl a y , a z ), fi is the chemical potential, and the 
interlayer coherent hopping has the form, 



which is strongly anisotropic and follows the theoret- 
ical predictions 2 ^. In particular, this momentum de- 
pendent form (2) has been experimentally verifie d 14 ' 15 . 
For this t-t'-J model (1), it has been argued that cru- 
cial requirement is to impose the electron single occu- 
pancy local constraint for a proper understanding of the 
physical properties of cuprate superconductors. To in- 
corporate the electron single occupancy local constraint, 
the charge-spin separation (CSS) fermion-spin theory has 
been proposed 2 ^, where the constrained electron opera- 
tors are decoupled as, C mT = h\ a ^S^ a and C iai = h\ a[ Sf a , 
with the spinful fermion operator hi aa — e~ l * iaa 7i ia rep- 
resents the charge degree of freedom together with some 
effects of the spin configuration rearrangements due to 
the presence of the doped hole itself (dressed holon), 
while the spin operator Si a represents the spin degree 
of freedom, then the bilayer t-t'-J Hamiltonian (1) can 
be expressed in this CSS fermion-spin representation as, 

H = t ^2(h i+ - a ^hi a iSfS i+fja + h\ +fjai h iai S ia S+ +fja ) 

i ff a 

- t' y^Xh+tar^i^^ta^i+fa + h \+Tal h ™l S ia S t+fa) 
if a 

i 

+ h \2i h inS a S+ + h\ n h i2i S^S+) - h \aa h ^a 

iaa 

+ J c ff Sia • Si+r)a + J c ff _L &il ' Si2, (3) 
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£_L.(k) = -i(cos k x — cos ky) 2 , 



(2) 



where J e g — J(l — S) 2 , J c ffj_ = «/±(l — S) 2 , and 8 — 
(h\ acr hi aa ) — (h\ a hi a ) is the doping concentration. It has 
been shown that the electron single occupancy local con- 
straint is satisfied in analytical calculations within this 
CSS fermion-spin theory, and the double spinful fermion 
occupancy are ruled out automatically 25 . Although in 
common sense hi aa is not a real spinful fermion, it be- 
haves like a spinful fermion 25 . As in the single layer 
case^, the kinetic energy terms in the bilayer t-t'-J 
model have been transferred as the dressed holon-spin 
interactions, which can induce the dressed holon pair- 
ing state (hence the electron Cooper pairing state) by 
exchanging spin excitations in the higher power of the 
doping concentration. Before calculation of the electron 
normal and anomalous Green's functions of the bilayer 
system in the SC state, we firstly introduce the SC order 
parameter. As we have mentioned above, there are two 
coupled Cu02 planes in the unit cell, and in this case, 
the SC order parameters for the electron Cooper pair is 
a matrix A = + <j x At, with the longitudinal and 
transverse SC order parameters in the CSS fermion-spin 
theory can be expressed as, 

Al = (c ia ^c i+fjai - cj al cj +f)a1 ) 

= {hia]h l+fj aiS+ a S l+fia - h ia ih l+fja] S la S+ +fia ) 

= -XiAfci, (4a) 

A T = (C| 1T Cj 2i - C iXi C m ) 

= {h l i^h i2 iSf l S~ 2 - hmhaiSaSe) 

= -x±&hT, (4b) 
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respectively, where the spin correlation functions Xi = 
( S ia S i~+f)a) and X± = (Sti^a), and the longitudinal 
and transverse dressed holon pairing order parameters 
= (hi+fiaihiai — h i+ f ia ]hia\) and A^t = (^i2jAiT — 
hi2-[hni). 

Within the t-J type model, robust indications 
of superconductivity with the d-wave symmetry in 
doped cuprates have been found by using numerical 
techniques 2 ^. On the other hand, it has been argued that 
the SC transition in doped cuprates is determined by the 
need to reduce the frustrated kinetic energy 2 ^. Although 
it is not necessary for the strong coupling of the electron 
quasiparticles and a pairing boson in their arguments 27 , 
a series of the inelastic neutron scattering experimental 
results provide a clear link between the electron quasi- 
particles and magnetic excitation s 28 ' 29 . In particular, an 
impurity-substitution effect on the low energy dynamics 
has been studied by virtue of the ARPES measurement^, 
this impurity-substitution effect is a magnetic analogue 
of the isotope effect used for the conventional super- 
conductors. These experimental results^ reveal that 
the impurity-induced changes in the electron self-energy 
show a good correspondence to those of the magnetic 
excitations, indicating the importance of the magnetic 
fluctuation to the electron pairing in cuprate supercon- 



ductors. Recently, wei 9 . have developed the kinetic en- 
ergy driven SC mechanism based on the CSS fermion- 
spin theory 2 ^, where the dressed holons interact occur- 
ring directly through the kinetic energy by exchanging 
spin excitations, leading to a net attractive force between 
dressed holons, then the electron Cooper pairs originat- 
ing from the dressed holon pairing state are due to the 
charge-spin recombination, and their condensation re- 
veals the SC ground-state. Within this SC mechanism 19 , 
the doping and temperature dependence of the electron 
spectral function of the single layer cuprate supercon- 
ductors in the SC state has been discussed^. In this 
section, our main goal is to generalize these analytical 
calculations from the single layer case to the bilayer sys- 
tem. As in the case for the SC order parameter, the full 
dressed holon normal and anomalous Green's functions 
can also be expressed as g(k, lj) = <?i,(k, u>) + a x gT(k, oS) 
and Syt(k,a;) = Q^(k,o;) + cr^^k, u>), respectively. We 
now can follow the previous discussions for the single 
layer cas o 18 ' 19 , and evaluate explicitly these correspond- 
ing longitudinal and transverse parts of the full dressed 
holon normal and anomalous Green's functions as [see 
the Appendix], 



( } " 2 J- hFA V^e^ + ^TE^rJ ' 

v— 1,2 v / 

^(M-^E^f- — \ 

LK ' ' 2 ^ hFA 2E hvk \cj-E hvk uj + E hvk )' 

y— 1,2 x ' 

4(k )W) = --v (-ir+^fMM ( i i ) , 



where the dressed holon quasiparticle coherence fac- 
tors U 2 huk = [1 + U/E hvk }/2 and = [1 - 

Cwk/-2'/u/k]/2, the dressed holon quasiparticle disper- 
sion EhuU = \J[fivk] 2 + | A^(k) | 2 , the renormalized 

dressed holon excitation spectrum £„k = Z hFA^^ with 
the mean-field (MF) dressed holon excitation spectrum 
6* = Zt X i7k - Zt' X2 y k -n+ (-If +1 X ±t±(k), where 
the spin correlation function \i — (SiaS i+fa ) , 7k = 
(l/^)E^e 4k ^, 7 'k = (V^)Ef e jk ' f , Z is the num- 
ber of the nearest neighbor or next nearest neighbor 
sites, the renormalized dressed holon pair gap func- 
tion A$(k) = Z$ A [A hL (k) + (-l)^A ftT (k)], with 
v = 1 ( v = 2) for the bonding (antibonding) case, 



where A M (k) = S^fk.w) | u=0 = K hL ^\ with 7 ^ = 
(cosfca, - cosk y )/2, A hT (k) = E^(k,a;) | w =o= A hT , the 
dressed holon quasiparticle coherent weights z\^p A = 

Z hFl _ Z hF2> Z hFA = Z hFl + Z hF2> Wltn Z hFl 

1 - S^ o) (k ,w) | w=0 , and Z^ 2 = E^ o) (k ,w) | w=0 , 

where ko = [ir, 0] , (k, u>) and E^ "* (k, u) are the cor- 
responding antisymmetric parts of the longitudinal and 
transverse dressed holon self-energy functions E^ (k, w) 
and E^j?(k, u>), while the longitudinal and transverse 
parts of the dressed holon self-energy functions E^ (k, oS) 
and E^ (k, lu) have been evaluated as, 
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-Rp+ q+k i ^ 5t(p + k, ip m + icj„)n TL (p, q, ip m )], 



(6a) 



7^2 H^P+q+k^ H 5^(P + k , iPm + iw„)II TT (p, q, ip m ) 



p,q 
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-Rp+ q+k i ^2 9l(p + k, ip TO + iw n )n iT (p, q, ip m )], 



(6b) 



Ilii(k,w n ) = ^^[i?W q+k I^Q[(_ p _k,-ip m -^„)n LL (p,q,ip m ) 



p.q 
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fip+ q+k -5 ^2 S t(~P ~ k, -ip m - iuj n )U TL (p, q, ip m )], 



72 Xl^P+q+k^ XI 3 t(-P - k , -ij> m - iw„)n TT (p, q, zp„ 



iV 2 
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p+q+k 



-iPn 



iw„)IlLT(p,q, ipr; 



(6c) 



(6d) 



where i?^ = [Z(i 7k - i' 7 ^)] 2 + ^(k), = 2Z(t 7k 
*'7k)*J-(k), and the spin bubbles H^rf (p, q, ip m ) 



(!//3) J2 lqm D n ] (q. km)D%> (q + p, jg TO + ip m ), with 



where = A(Ai 7k -A 2 ) -A' (2xl7'k-X2) - J e ffj_[X-L + 
2xl(-l)1[e±(k) + (-l)1, Ai = 2e|| X f+Xi, A 2 = e ,, Xl + 
2 X f, A = 2ZJ cff , X' = AZfot', en = l + 2t^/J eS , e±(k) = 
L, T and V = L, T, and the MF spin Green's function } + ^^txCk)/^ the spin correlation functions x f = 
i?W(k, W ) = ^(k.c) + a x 4°)(k, W ), with the cor- * 2 = xl = the dressed 

responding longitudinal and transverse parts have been holon particle- hole order parameters (f>i — (h^hi+fiaa) , 
given byS 2 ., <t>2 = {h\ aa h i+ f ar7 ) , <f>x = {h\ la h i2a ) , and the MF spin 

excitation spectrum, 



4 0) (k,c) = \ £ 



B. 



i/k 



2 ^ w 2 -c^ 2 ' 

y=l,2 



i>=1,2 ^ k 



(7a) 
(7b) 



r , 2 - X 2 



M - a£||Xi7k - ^ ae ||Xi) (1 - £||7k) + M 3 - -|axi - c*Xi7k 



-7k) 



+ A 



/2 



" i X27 k 22~ X2 J 7 k + 2 V z" X2 



AA'c 



Xi(l ~ £||7k)7'k + ^(XlVk ~ C2)(e|| - 7k) 



1 



+ 7k(C|-e t |X^7 k )--e||(C 2 -X27k) 



AJ cff± a ^ -e x (k)(e|| - 7k )[Cj_ + Xi(-l)1 



+ (1 - e ,| 7k )[Ci + X fex(k)(-l)1 + [e x (k) + (-1)"] 



2 e ||(Cx-x±7k) + (Ci-e||xi7k)(-ir 



+ A'J off± a ^ 7 ' k [C"l + x!ex(k)(-l)1 - -e±(k)[C' ± + x 2 (-l)"] 



2(x±7'k-^)+xl7'k(-l)' y 



[e ± (k) + (-1)"] [> + i J c 2 ffJ .Mk) + (-If] 2 , 



(8) 



where A 3 = aCi + (1 -a)/2Z, A 4 = aCf + (1 -a)/4Z, tion functions Ci = 0-/ Z 2 ) T,irf> (S+^S' J , 
A$ = aC3 + (1 — a)/2Z, and the spin correla- 
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C 2 - (l/^ 2 ) J2fjt ( S i+f,a S i+Ta)> ^3 

W&m^si^S'^), c ± = (i/z)z f1 (s+sr +fi2 }, 

and C'l = (l/Z)J2 f ( S a S Ut2)- In order to satisf y the 
sum rule of the spin correlation function (S^S^) = 1/2 
in the case without AFLRO, the important decoupling 
parameter a has been introduced in the above calcu- 
lation as in the single layer cas o 18 i 19 ' 22 , which can be 
regarded as the vertex correction. 

With the help of the longitudinal and transverse parts 



Gi(k, lu) = 



+ 

G T (k,w) = 



+ 

lt(k,a;) = 



4(k,w) = 



where L$(k,p) = [coth(/?w„ p /2) - th(/3£, w _ k /2)]/2 

and L$(k,p) = [cothO^p/2) + th(/3£ h( , p _ k /2)]/2, 
then the longitudinal and transverse parts of the 
electron spectral function Aj,(k, 10) = — 2ImGL(k, us) 



of the full dressed holon normal and anomalous Green's 
functions in Eq. (5) and MF spin Green's function 
in Eq. (7), we now can calculate the electron nor- 
mal and anomalous Green's functions G(i — j,t — t') = 
((C la (t);C%(t'))} = G L (i-j, t-t') + a x G T (t -j,t- f ) 
andrt(i-j,t-t') = ((c} r (t);C\ l (t')))=Tl(i-j,t-t') + 

a x Tj,(i—j, t—t'), where these longitudinal and transverse 
parts are the convolutions of the corresponding longitudi- 
nal and transverse parts of the full dressed holon normal 
and anomalous Green's functions and MF spin Green's 
function in the CSS fermion-spin theory, and can be eval- 
uated explicitly as, 



(9a) 



(9b) 



(9c) 



(9d) 

I 

and Ay(k, w) = — 2ImGT(k, uj) and SC gap func- 
tion A L (k) = (l//?)E l c„r[(k,^„) and A T (k) = 
(V/3) J2iu) n r^(k, iui n ) are obtained as, 
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A^(p-k)^ r 4i)(kip) ^ 



v-v- 7 (m) A^(p-k) B, 
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1 



'hfip— k t^vp 



v w — Eh^p-k + u> U p 



w + Ehup-u. — u V p ) ' ^ \k> — Eh^ji-k — u) V p w + Ehup-u. + u) V p 

EE(-D^< ) / 2 | (p " k) g^ k } (k, P ) ( - g 1 +lj 

p ~J ^h/ip-k ^vp |_ \W — ^hup-k + W„ p 

— — ) + L«(k, P ) ( — p 1 - — Y 

W + &h/ip-k — ^up J \W — -C^p-k — Wj/p W + &h/ip-k + U V p) _ 



1 

8iV 



I 

Ai(k, ^) = 41 — {^( k . P)[^ P -k<5(^ + *Wk - ^p) + ^p- k ^(^ - ^VP-k + ^ P )] 

P 

+ ^( k »P)[^ MP -k^(^ + ^VP-k + ^ P ) + V£ w _ k % - Sfc/.p-k - w«p)]}, (10a) 
A T (k, W ) = vr-l-E EC- 1 )^"^^^^' P)[^P-k^( w + ^VP-k - ^ P ) + V^ p _ v *( W - £ VP -k + ^ P )] 

P 

+ i| I i ) (k,p)[C/ 2 M p_ k ^(^ + ^vp-k + u vv ) + V^ p _ k 5(u - E hMP _ k - ^ p )]}, (10b) 
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16A ^ hFA 
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th[-/3£ Vp _ k ]coth[-/3o 



VpJ 



(10c) 
(lOd) 



With the above longitudinal and transverse parts of the 
SC gap functions in Eqs. (10c) and (lOd), the corre- 
sponding longitudinal and transverse SC gap parameters 
are obtained as A^ = ~xi^hL and At = —Xx^hT, 
respectively. In the bilayer coupling case, the more ap- 
propriate classification is in terms of the spectral func- 
tion and SC gap function within the basis of the an- 
tibonding and bonding component o 13 ' 14 i 15 ' 16 i 17 . In this 
case, the electron spectral function and SC gap parame- 
ter can be transformed from the plane representation to 
the antibonding-bonding representation as, 



A (a \k,u) 

A (6) (k,w) 

A (a) 
A (b) 



-A t Ql,u>) 
+ A T (k,io) 



A T , 
At- 



(11a) 

(lib) 

(11c) 
(lid) 



respectively, then the antibonding and bonding parts 
have odd and even symmetries, respectively. 



III. ELECTRON STRUCTURE OF BILAYER 
CUPRATE SUPERCONDUCTORS 

We now begin to discuss the effect of the bilayer in- 
teraction on the electronic structure in the SC state. We 
first plot, in Fig. 1, the antibonding (solid line) and 
bonding (dashed line) electron spectral functions in the 
[ir, 0] point for parameters t/J = 2.5, t'/t — 0.3, and 
t±/t= 0.35 with temperature T — 0.002 J at the doping 
concentration 8 — 0.15. In comparison with the single 
layer case^, the electron spectrum of the bilayer system 
has been split into the bonding and antibonging compo- 
nents, with the bonding and antibonding SC quasipar- 
ticle peaks in the [ir, 0] point are located at the differ- 
ent positions. In this sense, the differentiation between 
the bonding and antibonding components of the electron 
spectral function is essential. The antibonding spectrum 
consists of a low energy antibonding peak, corresponding 
to the SC peak, and the bonding spectrum has a higher 
energy bonding peak, corresponding to the hump, while 
the spectral dip is in between them, then the total con- 
tributions for the electron spectrum from both antibond- 
ing and bonding components give rise to the PDH struc- 
ture. Although the simple bilayer t-t'-J model (1) can- 
not be regarded as a comprehensive model for a quanti- 
tative comparison with bilayer cuprate superconductors, 
our present results for the SC state are in qualitative 
agreement with the major experimental observations on 
bilayer cuprate super conduct or 8,9, 10,11, i6 _ 

We now turn to discuss the doping evolution of the 
electron spectrum of bilayer cuprate superconductors in 
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FIG. 1: The antibonding (solid line) and bonding (dashed 
line) electron spectral functions in the [tt, 0] point for t/J — 
2.5, t'/t = 0.3, and t±/t = 0.35 with T = 0.002J at S = 0.15. 
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FIG. 2: The electron spectral functions at [n, 0] point for 
t/J = 2.5, t'/t = 0.3, and t±/t = 0.35 with T = 0.002J at 
S = 0.09 (solid line), 5 = 0.12 (dashed line), and 5 = 0.15 
(dotted line). 



the SC state. We have calculated the electron spec- 
trum at different doping concentrations, and the result 
of the electron spectral functions in the [tt, 0] point for 
t/J = 2.5, t'/t = 0.3, and t x /t = 0.35 with T = 0.002J 
at 5 — 0.09 (solid line), S = 0.12 (dashed line), and 
6 = 0.15 (dotted line) are plotted in Fig. 2. In compari- 
son with the corresponding ARPES experimental results 
of the bilayer cuprate superconductor Bi2Sr2CaCu20s+5 
in the SC state in Ref.— , it is obviously that the doping 
evolution of the spectral weight of the bilayer supercon- 
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FIG. 4: The antibonding (solid line) and bonding (dashed 
line) gap parameters as a function of the doping concentration 
with T = 0.002J for t/J = 2.5, t'/t = 0.3, and t±/t = 0.35. 



ductor Bi2Sr2CaCu208+,5 is reproduced. With increas- 
ing the doping concentration, both SC peak and hump 
become sharper, and then the spectral weights increase 
in intensity. Furthermore, we have also calculated the 
electron spectrum with different temperatures, and the 
results show that the spectral weights of both SC peak 
and hump are suppressed with increasing temperatures. 
Our these results are also qualitatively consistent with 
the ARPES experimental results on bilayer cuprate su- 
perconductors in the SC state^^. 

To better perceive the anomalous form of the antibond- 
ing and bonding electron spectral functions as a function 
of energy to for k in the vicinity of the [tt, 0] point, we 
have made a series of calculations for the electron spec- 
tral function at different momenta, and the results show 
that the sharp SC peak from the electron antibonding 
spectral function and hump from the bonding spectral 
function persist in a very large momentum space region 
around the [tt, 0] point. To show this point clearly, we 
plot the positions of the antibonding peak and bonding 
hump in the electron spectrum as a function of momen- 
tum along the direction [— 0.2tt,tt] — » [0,7r] — > [0.27r,7r] 
with T = 0.002J at 8 = 0.15 for t/J = 2.5, t'/t = 0.3, 
and t±/t = 0.35 in Fig. 3. Our result shows that there 
are two branches in the quasiparticle dispersion, with 
upper branch corresponding to the antibonding quasi- 
particle dispersion, and lower branch corresponding to 
the bonding quasiparticle dispersion. Furthermore, the 
BS reaches its maximum at the [V, 0] point. Our present 
result also shows that in analogy to the two flat bands ap- 
peared in the normal stated, both electron antibonding 
peak and bonding hump have a weak dispersion around 
the [tt, 0] point, in qualitative agreement with the ARPES 
experimental measurements on bilayer cuprate supercon- 
ductors in the SC state2££&±&±ii±l 

In the above calculations, we find that although the 
antibonding SC peak and bonding hump have different 
dispersions, the transverse part of the SC gap param- 
eter At ~ 0. To show this point clearly, we plot the 
antibonding and bonding gap parameters in Eqs. (11c) 



and (lid) as a function of the doping concentration with 
T = 0.002J for t/J = 2.5, t'/t = 0.3, and t±/t = 0.35 in 
Fig. 4. As seen from Fig. 4, both antibonding and bond- 
ing gap parameters have the same d-wave SC gap mag- 
nitude in a given doping concentration, i.e., A a rs A b . 
This result shows that although there is a single elec- 
tron interlayer coherent hopping (2) in bilayer cuprate 
superconductors in the SC state, the electron interlayer 
pairing interaction vanishes. This reflects that in the 
present kinetic energy driven SC mechanism, the weak 
dressed holon-spin interaction due to the interlayer co- 
herent hopping (2) from the kinetic energy terms in Eq. 
(3) does not induce the dressed holon interlayer pair- 
ing state by exchanging spin excitations in the higher 
power of the doping concentration. This is different from 
the dressed holon-spin interaction due to the intralayer 
hopping from the kinetic energy terms in Eq. (3), it 
can induce superconductivity by exchanging spin excita- 
tions in the higher power of the doping concentration^. 
Our this result is also consistent with the ARPES ex- 
perimental results of the bilayer cuprate superconductor 
Bi(Pb)2Sr2CaCu208+ J 14 ' 16 , where the SC gap separately 
for the bonding and antibonding bands has been mea- 
sured, and it is found that both d-wave SC gaps from 
the antibonding and bonding components are identical 
within the experimental uncertainties. 

To our present understanding, two main reasons why 
the electronic structure of bilayer cuprate superconduc- 
tors in the SC state can be described qualitatively in 
the framework of the kinetic energy driven supercon- 
ductivity by considering the bilayer interaction are as 
follows. Firstly, the bilayer interaction causes the BS, 
this leads to that the full electron normal (anomalous) 
Green's function is divided into the longitudinal and 
transverse parts, respectively, then the bonding and an- 
tibonding electron spectral functions (SC gap functions) 
are obtained from these longitudinal and transverse parts 
of the electron normal (anomalous) Green's function, 
respectively. Although the transverse part of the SC 
gap parameter At ~ 0, the antibonding peak around 



8 



the [tt, 0] point is always at lower binding energy than 
the bonding peak (hump) due to the BS. In this sense, 
the PDH structure in the bilayer cuprate superconduc- 
tors in the SC state is mainly caused by the BS. Sec- 
ondly, the SC state in the kinetic energy driven SC 
mechanism is the conventional BCS like as in the sin- 
gle layer cas o 18 ' 19 . This can be understood from the 
electron normal and anomalous Green's functions in Eq. 
(9). Since the spins center around the [ir, tt] point in 
the MF leve l 18 ' 19 ' 22 , then the main contributions for the 
spins comes from the [tt, tt] point . In this case, the 
longitudinal and transverse parts of the electron nor- 
mal and anomalous Green's functions in Eq. (9) can 
be approximately reduced in terms of u> vp= [„^ ~ and 
one of the self-consistent equations^ 2 - 1/2 = (5*^5" ) = 
V(4A0 E,,k(^k/^k)coth[(l/2)/3^ k ] as, 



G L (k,w)«i Z 



2 FA \ lj - £U ' w + E u]i 



ir, 



(12a) 



i/=l,2 



v. 



//k 



(12b) 



rt(k, w ) = - x: z 



1 ^ ^)A^(k) / 1 



2 ^ 9 FA ~ 2E„ k \u-E vi[ u + E vV 

v— 1,2 x 



4(k ) .) = i^(-i ) ^4i^ 



j/=l,2 



w - £7„k w + E,, 



(12c) 



(12d) 



where the electron coherent weig hts Z { p\ = Z$ A /2, the 
electron quasiparticle coherence factors /7 2 k ~ Vfi vk _ k 



and ^ 2 k 



U 2 



k , the SC gap function Aa (k) « 

hz ^ — k^) and the electron quasiparticle spectrum 
£vk ~ £^hi/k-k A , with kA = [7r,7r]- As in the sin- 
gle layer cas o 18 ' 19 , this reflects that the hole-like dressed 
holon quasiparticle coherence factors Vhuk and Uhvk and 
hole-like dressed holon quasiparticle spectrum Eh v \t have 
been transferred into the electron quasiparticle coher- 
ence factors t/^k and V uk and electron quasiparticle spec- 
trum £vk, respectively, by the convolutions of the corre- 
sponding longitudinal and transverse parts of the MF 
spin Green's function and full dressed holon normal 
and anomalous Green's functions due to the charge-spin 
recombination 27 . As a result, these electron normal and 
anomalous Green's functions in Eq. (12) are typical bi- 
layer BCS like 6 .. This also reflects that as in the single 
layer cas o 18 ' 19 , the dressed holon pairs condense with the 
d-wave symmetry in a wide range of the doping concen- 
tration, then the electron Cooper pairs originating from 
the dressed holon pairing state are due to the charge- 
spin recombination, and their condensation automati- 
cally gives the electron quasiparticle character. These 
are why the basic bilayer BCS formalism 6 is still valid in 
discussions of SC coherence of the quasiparticle peak and 
hump, although the pairing mechanism is driven by the 



intralayer kinetic energy by exchanging spin excitations, 
and other exotic magnetic scattering 28-29 is beyond the 
BCS formalism. 



IV. SUMMARY AND DISCUSSIONS 

We have studied the electronic structure of bilayer 
cuprate superconductors in the SC state based on the 
kinetic energy driven SC mechanism^ 9 -. Our results show 
that the electron spectrum of bilayer cuprate supercon- 
ductors is split into the bonding and antibonding com- 
ponents by the BS, then the observed PDH structure 
around the [tt, 0] point is mainly caused by this BS, with 
the SC peak being related to the antibonding compo- 
nent, and the hump being formed by the bonding com- 
ponent. The spectral weight increases with increasing 
the doping concentration. In analogy to the two flat 
bands appeared in the normal state, the antibonding 
and bonding quasiparticles around the [tt, 0] point dis- 
perse weakly with momentum, in qualitative agreement 
with the experimental observation on the bilayer cuprate 
superconductors^ 7 -' 8 -' 9 - ' 10 ' 11 . Our these results also show 
that the bilayer interaction has significant contributions 
to the electronic structure of bilayer cuprate supercon- 
ductors in the SC state. 

It has been shown from the ARPES experiments 2 -^ 
that the BS has been detected in both normal and SC 
states, and then the electron spectral functions display 
the double-peak structure in the normal state and PDH 
structure in the SC state. Recently, we 2 ^ have studied the 
electron spectrum of bilayer cuprate superconductors in 
the normal state, and shown that the double-peak struc- 
ture in the electron spectrum in the normal state is dom- 
inated by the BS. On the other hand, although the anti- 
bonding and bonding SC peaks have different dispersions, 
the antibonding and bonding parts have the same d-wave 
SC gap amplitude as mentioned above. Incorporating our 
previous discussions for the normal state case^ 2 - and the 
present studies for the SC state case, we therefore find 
that the one of the important roles of the interlayer co- 
herent hopping (2) is to split the electron spectrum of the 
bilayer system into the bonding and antibonding compo- 
nents in both normal and SC states. As a consequence, 
the well pronounced PDH structure of bilayer cuprate su- 
perconductors in the SC state and double-peak structure 
in the normal state are mainly caused by the BS. 
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APPENDIX A: DRESSED HOLON BCS TYPE 
NORMAL AND ANOMALOUS GREEN'S 
FUNCTIONS IN BILAYER CUPRATE 
SUPERCONDUCTORS 

In the single layer case, it has been showni^ that 
the dressed holon-spin interactions from the kinetic en- 
ergy terms of the t-t'-J model are quite strong, and 
in the case without AFLRO, these interactions can in- 
duce the dressed holon pairing state (then the electron 
Cooper pairing state) by exchanging spin excitations in 
the higher power of the doping concentration. Following 
their discussion s 18 ! 19 , we obtain in terms of Eliashberg's 
strong coupling theory^ that the self-consistent equa- 
tions that satisfied by the full dressed holon normal and 
anomalous Green's functions in the bilayer system in the 
SC state as, 



^(k.w) 



gWQL,w) + gWQi,w) [s < h > (k, w ) ff (k, w) 

- E^ ) (-k,-w)9 st (k,w)] , (Ala) 
0<°>(-k, -lu) M h \-K -u/)3t(_k, - U ) 



-k,-w) ff (k,w) , (Alb) 



respectively, where the MF dressed holon normal Green's 
function 2 ? 5 (°)(k,w) = s[ 0) (k,w) + a^Qs, w), with 
the longitudinal and transverse parts are evaluated as 
5 [ 0) (k, w ) - (l/2)E,=i,2(^-6k)^ and^(k,w) = 
(1/2) E^x^- 1 )^ 1 ^ - 6k)" 1 , respectively, while 
the dressed holon self-energy functions E^ (k, lu) — 
E^(k,w) + a x E$Qt,w) and E£ fc) (k,w) = E^(k,w) + 
<r x E^(k, a;), with the corresponding longitudinal and 
transverse parts have been given in Eq. (6). 

In the previous discussions of the electronic struc- 



ture for the single layer cuprate superconductors in the 
SC state 18 , it has been shown the self-energy function 

Ej (k, lu) describes the effective dressed holon pair gap 
function, while the self-energy function E^ (k, lu) de- 
scribes the quasiparticle coherence. Since E^ 1 ' (k, lu) is 
an even function of u, while E^(k, lu) is not, therefore 
for the convenience, the self-energy function E^ (k, lu) 
can be broken up into its symmetric and antisymmetric 
parts as, E<' l) (k,w) = E^ ) (k,cj)+wE^ ) (k,cj), then both 

E[g (k, lu) and E^(k, lu) are even functions of lu. Now 
we can define the dressed holon quasiparticle coherent 
weights in the present bilayer system as Z^ F1 (k, lu) — 

1 - Y,[ h L o) (k,Lo) and Z^,(k,w) = E^ o) (k,w). As in the 
single layer case^, we only discuss the low-energy behav- 
ior of the electronic structure of bilayer cuprate super- 
conductors, which means that the effective dressed holon 
pair gap functions and quasiparticle coherent weights 
can be discussed in the static limit, i.e., A^(k) = 

s£°(k,w) \ u=0 = Afc£(k) + a x A hT (k), Z-^(k) = 1 - 

E^ o) (k,a,) \ u=0 and Z^(k) = E^ o) (k, W ) | w=0 . As in 
the single layer case^, although ZhFiO*) and ZhFi (k) 
still are a function of k, the wave vector dependence 
may be unimportant. This followed from the ARPES 
experiments- that in the SC-state of bilayer cuprate su- 
perconductors, the lowest energy states are located at 
the [it, 0] point, which indicates that the majority con- 
tribution for the electron spectrum comes from the [tt, 0] 
point. In this case, the wave vector k in Z^fi0<) and 
ZhF2^) can be chosen as Z^p X = 1— E^^ (k) |k=[>r,o] an d 

Zfop 2 = (k) |k=[7r,o]- With the help of the above dis- 
cussions, the corresponding longitudinal and transverse 
parts of the dressed holon normal and anomalous Green's 
functions in Eqs. (Ala) and (Alb) now can be obtained 
explicitly as, 



n fir , a — 1 ( ^hvk -x- Yhnk \ 

9LiKUj) - 2 J- ZhFA \^E^ + ^Te^J ' 

^ ~[ 2 \LU — &huk LU + &hvk J 

^(k, W ) = --Y z^ ^M ( 1 _ 1 ) , 

2 ~( 2 ZEhvk \lu — Ehvk ^ + Ehuk ) 

S>t fk lu) - 1 V ( lY^Z^ A ^ (k) ( 1 1 ^ 

i/— 1,2 x ' 



with the dressed holon effective gap parameters and equations, 
quasiparticle coherent weights satisfy the following four 
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AhL = - 



Z 
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hFA 



7(2) 
J hFA 



E E 



32N 3 ^ ^ 7k -p+q 

k.q.p v,v' ,v" 



C w »(k + q) 



ry{v") D R 



F^(q:P)+^l»(k,q,p) , F(l„(q,p) + ^ 4 J, y „(k,q,p) 



//qj 



p2 

-^/liVk 



A^ T = - 



1 



+ 



32A 3 



E E (-ir +i/+,y " +1 axk+q) 



7( w ") D D , „. 

Z/ hFA- p "'p- p ' / q a (" ) 



k.q.p iaz/' ,v" 



UJ 



FiH„ («!, P) + (k, q, P) , (q p) + F%, v „ (k, q p 



,(2) 



(3) 



,(4) 



(k) 



1 



E 2 



kV 



E 2 



V 2 E E 

q.p v,v>' ,v>" 



[1 + (-l)^'+^ +1 ]a^(p + k )-^ i3,y ' pi3,yq 



-V'pt-t-Vq 



#£l»(q,p) 



^ 2 X"(q>p) 



^[/"p-q+ko] 2 [ u 



-ETii/'p-q+ko^ 



+ 



= 1 



^l-(qp) 



^l»(q,p) 



•'i/q 



y"p-q+k J 



^'p 



■^q 



p-q+koJ 



32A 2 



E E t 1 - (-1)^ 



7("") R , R 

"^^(p + ko) ,lFA " p " q 



q.p v,v' ,v" 

^l»(qp) 



-V'pC-i-Vq 



(2) 



-(q,p) 



p 



J v<\ 



Ehv' 



p— q+koJ 



^U(qp) 



^tl"(q>p) \ 



i/'p 



//q 



+ ^/ lI y"p-q+k ] 2 [<< 



■^z-'q 



p-q+koJ 



(A3a) 



(A3b) 



(A3c) 



(A3d) 



where C vv „{k) = [Z(t lk - ty k ) + (-l^'i^k)] 2 , 
^iL»(q>P) = «s(^q) + «B(^ 1 y'p) + 2ri B (cj, yq )n s (a; 1/ /p), 
^l-(k.q.P) = [2n F (^ k ) - l][uv p - 



^l«(q,p) = 



"sf^Vp^/iW'k 
1 + n_B(wiyq) + raB(ciVp) + 2nfl(o; I /q)TiB(w v /p), 

^l»(k,q,p) = [2n F {E hv ,, k ) - l][uv p + c^ q ][l + 

nsC^vq) + riB^p)]/^, Hili^(q,p) = 
nF(£^"p^q+k )[nBK' P ) - ns^q)] + ns(w„q)[l + 

"flk'p)], ^V»(q,Pj = «F(-E^''p-q+ko)[«B(^q) ~ 



(3) 

n B (uv' P )] + n B (uv' P )[l + n B (ui ucl )}, HlJ, u „(q,p) = 
[1 - nF(-Ewp-q+ko)][l + fi B (w„ q ) + n B (w„' p )] + 

IiJjKq)liJjK- p ), -ff^V"(q-P) = "'F(-E ; ft ! /"p-q+ko)[l + 

n B (^uq)+n B (iJ u , p )]+nB(uJuq)nB(^u'p), and k = [7r,0]. 
These four equations must be solved self-consistently 
in combination with other equations as in the single 
layer cas e 18 ' 19 , then all order parameters, decoupling 
parameter a, and chemical potential /1 are determined 
by the self-consistent calculation. 
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